Deuteron-Nucleus Collisions in a Multi-Phase Transport Model 



Zi-wei Lin 

Physics Department, The Ohio State University, Columbus, Ohio 43210 



CO 

o 

O 

CN 

< 

OS 



> 
in 

CN 

o 

o 

co 
o 



o 

=5 



X 



Che Ming Ko 

Cyclotron Institute and Physics Department, Texas A&M University, College Station, Texas 77843 

(Dated: February 8, 2008) 

Using a multi-phase transport (AMPT) model, we study pseudo-rapidity distributions and trans- 
verse momentum spectra in deuteron-gold collisions at RHIC. We find that final-state partonic and 
hadronic interactions affect the transverse momentum spectrum of protons more than those of kaons 
or pions. Relative to p+p collisions at same center-of-mass energy per nucleon pair, the effect of 
final-state interactions on the charged particle transverse momentum spectra in d+Au collisions is 
much smaller than observed in experimental data, indicating that initial-state effects such as the 
Cronin effect are important. 

PACS numbers: 25.75.-q, 24.10.Lx 



I. INTRODUCTION 

To study heavy ion collisions at the Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven National Laboratory, 
in which a deconfined plasma of quarks and gluons is ex- 
pected to be formed, we have developed a multi-phase 
transport (AMPT) model that includes both final-state 
partonic and hadronic interactions Q, 0, ■ The model 
has been very useful for understanding various observ- 
ablcs in Au+Au collisions at RHIC such as the rapidity 
and transverse momentum distributions of various parti- 
cles pi 0, 0] as well as charmonium Q and strangeness 
[f| productions. In particular, it allows one to study both 
thermal and chemical equilibration in the partonic and 
hadronic matter formed in these collisions. The AMPT 
model has also been extended to include the string melt- 
ing mechanism 0, Q , in which soft strings produced from 
initial nuclcon-nucleon interactions are converted directly 
to partons, in order to explain the measured elliptic flow 
@ and two-pion correlation functions at RHIC. 

In spite of its success, the AMPT model has many un- 
certainties in its input physical parameters, particularly 
those related to the initial conditions introduced to the 
model. For example, the number of initial minijet par- 
tons in the AMPT model, that is given by the hard pro- 
cesses from the HIJING model 0, 0, p2j , depends on 
the nuclear shadowing, i.e., the modification of the par- 
ton distributions in a nucleon when it is in a nucleus. As 
a result, the final particle multiplicity produced in rela- 
tivistic heavy ion collisions is affected by nuclear shadow- 
ing 0, Q . The production of high particles or hadrons 
made of heavy quarks, which is described by perturba- 
tive QCD processes, is even more sensitive to the nuclear 
shadowing effect. Furthermore, other nuclear effects such 
as the Cronin effect , which set in already in p-A col- 
lisions, need to be included. The uncertainties in the ini- 
tial condition also exist in other transport models with 
partonic degrees of freedom 0, 0] as well as theoretical 
models such as the hydrodynamical mo del [l6Lfl7Lfll .lT^| 
and the QCD saturation model [H EJ iMlal^P 



Since final-state interactions (FSI) are expected to be 
less important in dcutcron-gold collisions than in colli- 
sions between heavy nuclei, both the Cronin effect and 
the nuclear shadowing effect (e.g., through dilepton mea- 
surements 25] ) can be better studied in these collisions. 
Improved knowledge on these effects are useful for mak- 
ing reliable theoretical interpretations of the observa- 
tions in heavy ion collisions at RHIC. In this paper, we 
use the AMPT model to study the global observables 
in deuteron-gold collisions such as the pseudo-rapidity 
distributions and the transverse momentum spectra. Be- 
cause of the small interaction volume in d+Au collisions, 
we use the default AMPT model, i.e., without the string 
melting mechanism, in the present study, as the initial 
energy density produced in these collisions is expected to 
be small. 

This paper is organized as follows. In Sec. EI we briefly 
review the AMPT model. Results from the AMPT model 
on deuteron-gold collisions at RHIC are then shown 
in Sec. IIIII with the charged particle pseudo-rapidity 
distributions and their centrality dependence given in 
Sec. IIII Al the transverse momentum spectra of different 
particles in Sec. IIII Bl and the effects due to modifica- 
tions of the string fragmentation parameters and nuclear 
shadowing in Sec. IIII CI A summary is then given in 
Sec. IIVI Finally, the effect of centrality selection on the 
centrality dependence of charged particle pseudo-rapidity 
distributions is discussed in the Appendix. 



II. THE AMPT MODEL 

The AMPT model is a hybrid model that consists of 
four components: the initial conditions, the parton cas- 
cade, the conversion from partonic to hadronic matter, 
and the hadron cascade. In the default AMPT model 
jHQEJ[J0l the initial conditions are generated from the 
HIJING model (version 1.383 for this study), 

which usually uses a Woods-Saxon radial shape for the 
colliding nuclei and introduces a parameterized nuclear 
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shadowing function that depends on the impact param- 
eter of the collision |icj . Interactions among minijct 
partons, which are produced from initial hard nucleon- 
nucleon interactions, are modeled by the Zhang's parton 
cascade (ZPC) ^ij. After partons stop interacting, they 
recombine with their parent strings, which are produced 
from initial soft nucleon-nucleon interactions, and frag- 
ment to hadrons according to the Lund string fragmenta- 
tion model |26l l27j . Dynamics of resulting hadronic mat- 
ter is then described by a hadronic cascade based on the 
relativistic transport (ART) model j2^. Final hadronic 
observables including contributions from the strong de- 
cay of resonances are determined when the hadronic mat- 
ter freezes out. 

In this study, we use the Hulthen wave function [2j| to 
model the structure of a deuteron [3(| : 

u(r) = Ce~ ar (1 - e-" r ) , (1) 

where r represents the relative distance between the 
proton and neutron in the deuteron, the normalization 
constant C is determined from L u(r) 2 dr = 1, and 
the small -D-wave contribution (~6%) to the deuteron 
wave function is neglected. With a — (4.38 fm) -1 and 
fi = (1.05 fm) -1 , the root-mean-square radius of the 
deuteron, yj (r 2 ) /2, is 2.0 fm, consistent with the mea- 
sured value. 
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FIG. 1: Pseudo-rapidity distributions of charged particles in 
d+Au collisions at ^/sjvjv = 200 GeV with centralities de- 
termined from Npart- Solid and dashed curves represent re- 
sults from the AMPT model and the HIJING model (without 
quenching) , respectively. 



III. AMPT RESULTS ON DEUTERON-GOLD 
COLLISIONS 

We have studied deuteron-gold collisions at ^/snn = 
200 GeV with both deuteron and gold having the same 
energy of 100 GeV per nucleon. In this study, "minimum 
bias" d+Au events are defined as those within the im- 
pact parameter range between and 12 fm, and they are 
separated into different centrality bins. 

A. Pseudo-rapidity distributions of charged 
particles 

Fig. ^ shows the pseudo-rapidity distributions of 
charged particles in the six centrality bins of 0-10%, 10- 
20%, 20-40%, 40-60%, 60-80%, and 80-100%. Here, the 
centrality is determined according to the value of N part , 
i.e., the total number of participants in both deuteron 
and gold nuclei, in each event. The pseudo-rapidity in 
this study is evaluated in the nucleon-nucleon center-of- 
mass frame, and negative rapidities correspond to the 
fragmentation region of the gold nucleus. Solid curves 
are results from the AMPT model, while dashed curves 
are those from the HIJING model without quenching. 
Note that the AMPT model without final-state partonic 
and hadronic interactions is equivalent to the HIJING 
model without jet quenching if the popcorn mechanism, 
controlled mainly by parameters MSTJ(12) and PARJ(5) 
in the PYTHIA/JETSET program 26], is treated in 



the same way. It is seen that the asymmetry in the 
pseudo-rapidity distributions, e.g., by comparing values 
of dN c h I drj at rj = —2 and r\ = +2, decreases as collisions 
become less central, and the pseudo-rapidity distribution 
for the most peripheral bin is almost symmetric. As ex- 
pected, final-state interactions in d+Au collisions have a 
smaller effect on charged particle pseudo-rapidity distri- 
butions than in the case of central heavy ion collisions at 
RHIC 0,0. 

We note that the appearance of small bumps around 
r) ~ —4 is mainly due to interactions of produced par- 
ticles with initial incoming nucleons in the gold nucleus, 
which has an initial half width of Rau/j in the longi- 
tudinal direction with Rau being the hard-sphere radius 
of the gold nucleus and 7 denoting the Lorentz boost 
factor of a nucleon in the nucleon-nucleon center-of-mass 
frame. If the initial half width is reduced by a factor of 5, 
the charged particle pseudo-rapidity distributions would 
be smooth around 77 ~ —4 and also have slightly lower 
values in the region of— 6 < r\ < —3. 

Fig. shows the pseudo-rapidity distributions of 
charged particles for "minimum bias" d+Au events. 
Solid curves are results from the AMPT model, dashed 
curves with open circles are AMPT results without final- 
state interactions, and dashed curves with filled circles 
are those without both final-state interactions and the 
popcorn mechanism for baryon-antibaryon production 
IE U> 13- We- see that both final-state interactions and 
the popcorn mechanism broaden the pseudo-rapidity dis- 
tributions, especially in the fragmentation region of the 
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FIG. 2: Pseudo-rapidity distributions of charged particles a) 
for "minimum bias" events when different interactions are in- 
cluded and b) for events from AMPT with iVp art = 2 and with 
N^ vt = 1. 



gold nucleus, leading thus to a moderate increase of the 
particle multiplicity in the region of —6 < r] < —3. 
We note that without final-state interactions and the 
popcorn mechanism, results from the AMPT model are 
equivalent to those from the HIJING model (without jet 
quenching) as same values of a and b in Eq. are used. 

One can separate events according to the number of 
participants from the deuteron, i.e., Np art . As an ex- 
ample, Fig. 03 ) shows the pseudo-rapidity distributions 
of charged particles from the AMPT model for events 
with one or two nucleon(s) from the deuteron that are 
involved in primary collisions. About 2/3 of "minimum 



bias" events from the AMPT model have A^ art = 2 with 
average values of N^ t at 10-4 and the mean impact pa- 
rameter at 4.7 fm. On the other hand, for events with 
1 the average value of N™. t is 2.8 and the mean 
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impact parameter is 7.8 fm. 

In Fig. |3 we show by the solid curve with circles the 
average impact parameter as a function of the average 
value of N part at each centrality bin. Here, the number 
of participants is defined as the number of initial (projec- 
tile and target) nucleons involved in primary collisions, 
i.e., not including those initial nucleons that interact with 
produced particles in the final state. The solid curve with 
no symbols in Fig. |3| represents Np alt . We see that N^ art 
increases toward the value of 2 rather quickly from pe- 
ripheral to central collisions. The dashed curve with open 
circles shows the number of participating nucleons that 
also include initial nucleons that are involved in final- 
state interactions. It is seen that final-state interactions 
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FIG. 3: The average value of impact parameter (in fm) and 
JVp art (solid curve) as a function of the total N part . See text 
for details. 



increase the number of participating nucleons by more 
than 40% for the most central events (0-10% centrality 
bin). However, since the energies involved in final-state 
interactions are much lower than those in the primary 
collisions, they contribute much less to particle produc- 
tion. In this study, we thus only include initial nucleons 
that are involved in primary collisions for determining 

Npart- 

The centrality dependence of particle multiplicity and 
transverse energy at mid-pseudo-rapidities (—0.5 < r\ < 
0.5) are shown Fig. 0] The solid curve and the 
dashed curve with open squares represent, respectively, 
dN c h/dr]/Np ai t and dE'r/dT]/~Np ax t, where £t includes 
the contribution from neutral particles. In this study, 
Et = EsmO with 8 being the polar angle, and E is de- 
fined as the kinetic energy for baryons, the total energy 
including the mass for anti-baryons, and the total energy 
for all other particles. We observe that both the par- 
ticle multiplicity and transverse energy per participant 
gradually decrease with increasing N par t- However, we 
shall see in Sec. lIII Cl and the Appendix that this moder- 
ate decrease with centrality will change when centralities 
are determined differently or when the string fragmenta- 
tion function is modified in more central d+Au collisions. 
Fig. 0| also shows that the ratio of the transverse energy 
and particle multiplicity per participant (solid curve with 
filled squares), i.e., the average transverse energy per 
charged particle, is rather flat. Also shown by the solid 
curve with triangles is the mean transverse momentum 
of charged particles, which is seen to change little with 
Npart as well. 

The centrality dependence of the charged particle mul- 
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FIG. 4: Centrality dependences of charged particle multiplic- 
ity and transverse energy (in GeV) per participant, the mean 
transverse energy (in GeV) per charged particle, and the mean 
Pt (in GeV/c) of charged particles at mid-pseudo-rapidity. 



(curve with filled diamonds). Furthermore, the decrease 
with centrality becomes stronger as the pseudo-rapidity 
changes from negative values to more positive values. 
These are consistent with the picture that multiple inter- 
actions in the gold nucleus push particle productions to- 
ward the negative rapidity region. Since we find that the 
mean transverse energy per charged particle at a given 
pseudo-rapidity does not change much with centrality, 
the centrality dependence of (d-Er/d^/Nparti i.e., the 
transverse energy per participant, has a similar behavior 
at each pseudo-rapidity as (dN c h/ drj) /N part . 



B. Transverse momentum spectra 

The transverse momentum spectra of charged parti- 
cles in the pseudo-rapidity range of — 1 < rj < 1 for both 
"minimum bias" and 0-20% central d+Au collisions are 
shown in Fig. |SJ Solid curves with circles are results from 
the AMPT model while dashed curves are those with- 
out final-state partonic and hadronic interactions. For 
comparisons, the AMPT results for minimum bias p+p 
collisions are shown by the solid curve. In all transverse 
momentum distributions and their ratios shown later, we 
include the statistical errors in the AMPT model, which 
become large at moderately high px- 
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FIG. 5: Centrality dependence of charged particle multiplicity 
per participant at different pseudo-rapidities. 



tiplicity per participant at different pseudo-rapidities is 
shown in Fig. [21 where the value at each r\ represents the 
average value within the pseudo-rapidity range of 77 + 0.5. 
We see an increase of this quantity at the backward 
pseudo-rapidity r/ = —4 (curve with open diamonds) and 
a fast decrease at the forward pseudo-rapidity rj = +4 



FIG. 6: Transverse momentum spectra of charged particles 
within — f < rj < 1 for d+Au collisions from AMPT with 
(solid curves with circles) or without (dashed) final-state in- 
teractions. The solid curve is the AMPT result for p+p col- 
lisions at the same center-of-mass energy per nucleon pair. 

To better examine the distributions shown in Fig. [HJ 
ratios of the pt spectra of charged particles from d+Au 
collisions over those from minimum bias p+p collisions 
are shown in Fig.[7| We see that, even without final-state 
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partonic and hadronic interactions, the ratios (dashed 
curves) are not flat but instead tend to increase with p^ 
in the range of — 3 GeV/c shown in the figure. This be- 
havior is mainly a result of different scalings for hard and 
soft processes with respect to the number of binary col- 
lisions in the initial conditions from the HIJING model. 
While particle production at high enough pt scales with 
the number of initial binary collisions, low momentum 
particles have a weaker dependence. Final-state inter- 
actions further modify the ratios of p^ spectra, as seen 
from comparisons between the solid and dashed curves. 

Recently, the transverse momentum spectra of charged 
particles in d+Au collisions have been measured at RHIC 
[31| . Compared with the experimental data at < 
Pt < 3 GeV/c, the pt dependence of the ratios from 
the AMPT calculations shown in Fig. [7J either with or 
without final-state interactions, are much weaker. Thus, 
interactions in the final state are not the main reason for 
the observed strong p^ dependence of the ratios of trans- 
verse momentum spectra, suggesting that initial-state ef- 
fects such as parton momentum broadening due to the 
Cronin effect [H H2 (not included in the AMPT model 
so far) are important in d+Au collisions. 



16 



V 

V 12 



E- 



T3 

O 



(no FSI, minimum bias)/pp 
(AMPT, minimum bias)/pp 
(no FSI, 0-20% central)/pp 
(AMPT, 0-20% central)/pp 




h+h" 



12 3 

p T (GeV/c) 

FIG. 7: Ratios of the charged particle transverse momentum 
spectra for d+Au collisions over the spectra for p+p colli- 
sions from AMPT with (solid) or without (dashed) final-state 
interactions. 

For identified particles such as ti + ,K + and protons, 
their transverse momentum spectra within the rapidity 
range of — 1 < y < 1 from "minimum bias" d+Au colli- 
sions are shown in Fig. [SJ with the proton spectra scaled 
down by a factor of 10. Solid and dashed curves rep- 
resent AMPT results with and without final-state inter- 
actions, respectively. Compared with results from the 
AMPT model on central heavy ion collisions, e.g., Fig. 2 
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FIG. 8: Transverse momentum spectra of tt + ,K + and pro- 
tons from AMPT with (solid) or without (dashed) final-state 
interactions for "minimum-bias" d+Au collisions. The proton 
spectra are scaled down by a factor of 10. 



of Ref. Q for central Pb+Pb collisions at SPS, the effect 
of final-state interactions on the transverse momentum 
spectra in d+Au collisions are much weaker. 

Fig. shows the ratios of the pt spectra from AMPT 
over those from the AMPT model without final-state in- 
teractions. Dashed curves represents the ratios of the 
"minimum bias" spectra shown in Fig. |SJ while solid 
curves represent the ratios of the spectra from 0-20% 
d+Au collisions. We see that final-state interactions 
modify the pt spectra of ir, K and protons differently. 
Proton spectra show the largest FSI effect with low Pt 
protons shifting to higher p^, and the ratio for 0-20% 
central events has a stronger j?t dependence than that 
for "minimum bias" events. Both features are consis- 
tent with the existence of some transverse flow in d+Au 
collisions. We find that the average number of partonic 
collisions per parton is about 0.06 (with 3 mb parton scat- 
tering cross section) and the average number of hadronic 
collisions per produced hadron is about 1 in "minimum 
bias" d+Au collisions from the AMPT model, compared 
with about 1 partonic collision per parton and 8 hadronic 
collisions per hadron in central Au+Au collisions at 200^4 
GeV. As a result, the transverse flow due to final-state 
interactions in d+Au collisions is much weaker than in 
central heavy ion collisions. We also note that in the 
default AMPT model used in this study, final-state in- 
teractions among partons are less important than those 
among hadrons as only minijet partons are included in 
the parton cascade. 

Ratios of the pt spectra in d+Au collisions over those 
in p+p collisions from the AMPT model are shown in 
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FIG. 9: Ratios of the transverse momentum spectra of it , 
and protons from AMPT over those from the AMPT 
model without final-state interactions for "minimum bias" 
(dashed) and 0-20% central (solid) d+Au collisions. 



FIG. 10: Ratios of the transverse momentum spectra from 
AMPT for d+Au collisions over those for p+p collisions. 
Solid and dashed curves correspond to 0-20% central and 
"minimum-bias" d+Au collisions, respectively. 



Fig. Similar to Fig. we observe the strongest pr 
dependence for protons. Also, ratios of the pt spectra of 
kaons have a stronger dependence than those of pions. 
We note that, since the initial-state parton broadening 
due to the Cronin effect has not yet been included 
in the AMPT model, the pt dependence of these ratios 
may get even stronger if the Cronin effect is taken into 
account. 



C. Effects of the string fragmentation parameters 
and nuclear shadowing 

In the default AMPT model, hadron production after 
the partonic phase is described by the Lund string frag- 
mentation model. In this model, the longitudinal mo- 
mentum distribution of a hadron with transverse mass 
tox produced from the string fragmentation is given by 
the following symmetric splitting function |27l |: 

f(z) oc z-\l - z) a exp(-6 m 2 Jz), (2) 

where z denotes the light-cone momentum fraction of the 
produced hadron with respect to that of the fragment- 
ing string. In the HIJING model, which reproduces the 
experimental charged particle multiplicities in high en- 
ergy p+p and pp collisions, the default values of a = 0.5 
and b = 0.9 GeV -2 are used. However, in order to re- 
produce the rapidity distributions of charged particles in 
central Pb+Pb collisions at the CERN-SPS energy using 
the AMPT model, we find that values of these param- 
eters need to change to a — 2.2 and b — 0.5 GeV -2 



|2j, |3(. The increase of a and decrease of b relative to 
the default values in the HIJING model soften the split- 
ting function and thus enhance the total charged particle 
multiplicity by about 20% in central Pb+Pb collisions at 
SPS. The change of these two parameters can perhaps 
be attributed to the modification of string fragmentation 
when multiple strings are produced and overlap in heavy 
ion collisions, although it is not known yet how these 
parameters should be modified in different colliding sys- 
tems. Since peripheral d+Au collisions are expected to 
behave similarly as p+p collisions, and even in central 
d+Au collisions the string density is much lower than in 
central heavy ion collisions, we have assumed that the 
parameters a and b of Eq. J5J for deuteron-gold collisions 
have the same values as in p+p collisions, i.e., they take 
the default values in the HIJING model. 

However, since the values of a and b parameters need to 
be modified in order to describe the total multiplicity in 
central heavy ion collisions, it is still possible that these 
parameters are also modified in non-peripheral deuteron- 
gold collisions. If in central deuteron-gold collisions the 
parameters a and b already have the modified values as 
in heavy ion collisions, the pseudo-rapidity distribution 
for the 20% most central deuteron-gold collisions (se- 
lected from Np ar t) would be given by the dashed curve 
in Fig. 1111 which shows a substantial increase compared 
with the previous AMPT result obtained with default a 
and b values in the HIJING model (solid curve). Com- 
parison of these results for deuteron-gold collisions with 
upcoming experimental results at RHIC will thus allow 
us to learn how the string fragmentation process is mod- 
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ified in different colliding systems or centralities. 
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FIG. 11: Charged particles pseudo-rapidity distributions for 
the most central 20% deuteron-gold collisions from AMPT 
with default values of a and b parameters (solid), with nuclear 
shadowing turned off (dot-dashed), or with modified a and b 
values as in heavy ion collisions (dashed). 
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FIG. 12: Ratios of the transverse momentum spectra of 
tt + ,K + and protons from AMPT with nuclear shadowing 
over those without nuclear shadowing for 0-20% central d+Au 
events. 



When the nuclear shadowing effect is turned off in the 
AMPT model, the pseudo-rapidity distribution for the 
20% most central deuteron-gold collisions is shown by the 
dot-dashed curve in Fig. 1111 which is not much different 
from the solid curve obtained with the nuclear shadow- 
ing. Thus, pseudo-rapidity distributions of charged par- 
ticles in deuteron-gold collisions are not very sensitive to 
the nuclear shadowing effect. This is different from that 
seen in collisions between heavy nuclei, where the effect 
is much larger as nuclear shadowing affects the produc- 
tion of minijet partons, which scales with the number of 
binary collisions. The effects of nuclear shadowing on pj> 
spectra of tt + ,K + and protons, shown in Fig. ^] for 0- 
20% central d+Au collisions, are also relatively small and 
do not have a strong dependence within < pt < 3 
GeV/c. Note that, although the nuclear shadowing effect 
is unimportant for global observables in small colliding 
systems such as deuteron-gold collisions, it affects sig- 
nificantly observables which are dominated by partonic 
interactions such as open charm production |25| or the 
yield of high px particles. 



IV. SUMMARY 

Using a multi-phase transport (AMPT) model that 
includes both final-state partonic and hadronic inter- 
actions, we have studied the pseudo-rapidity distribu- 
tions of charged particles, their centrality dependence, 



and transverse momentum spectra of different particles 
in deuteron-gold collisions at ^snn = 200 GeV. Due to 
the asymmetry of these collisions, the centrality depen- 
dence of charged particle multiplicity per participant is 
very different at different pseudo-rapidities, and it goes 
from increasing with centrality at the backward pseudo- 
rapidity region (the fragmentation region of the gold nu- 
clei) to decreasing with centrality at the forward region. 
The charged particle pseudo-rapidity distribution in cen- 
tral deuteron-gold collisions is also sensitive to the val- 
ues of the parameters used in the string fragmentation 
function. Furthermore, we find that, although final-state 
partonic and hadronic interactions modify the transverse 
momentum spectra of charged particles in d+Au colli- 
sions relative to scaled p+p collisions, the pt dependence 
of the modification due to final-state interactions is much 
weaker than observed at RHIC. Thus, initial-state effects 
such as parton momentum broadening due to the Cronin 
effect are important in deuteron-gold collisions. How- 
ever, we find that final-state interactions have a much 
stronger effect on the pt spectra of protons than those 
of kaons or pions. Comparison of these predictions with 
the experimental data will thus help us to learn more 
about initial-state effects on transverse momentum spec- 
tra and to study whether string fragmentation is mod- 
ified in deuteron-gold collisions as in central heavy ion 
collisions. 
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APPENDIX: EFFECTS OF DIFFERENT 
CENTRALITY SELECTIONS 

To extract useful information from the centrality de- 
pendence of global observables, it is helpful to study 
how the criteria of centrality selection affects the cen- 
trality dependence. The number of participants and the 
impact parameter are known in theoretical models, but 
they cannot be directly measured in experiments. In- 
stead, experimental centrality selections usually involve 
cutting on the charged particle multiplicity or the trans- 
verse energy. Thus, in addition to separating events into 
different centrality bins using the total number of par- 
ticipants, we have also tried three other methods, using 
the impact parameter, the number of charged particles 
within — 1 < rj < 1, and the total number of charged par- 
ticles. For example, when the total number of charged 
particles is used for the centrality selection, we order all 
35,000 events from the AMPT model by the total num- 
ber of charged particles in each event. The centrality bin 
of, say 0-10%, then consists of the first 3,500 events in 
that ordered list. 
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FIG. 13: The average value of impact parameter as a function 
of Np ar t with different criteria for the centrality selection, to- 
gether with the RMS width of the impact parameter in each 
centrality bin for two of the curves. 

In Fig. H3 we plot the average value of impact parame- 
ters versus the average value of N part at six centrality bins 
of - 10 - 20 - 40 - 60 - 80 - 100%, that are selected using 
these four different methods. The solid curve with circles 
corresponds to the centrality selection from N part , and it 
thus has the largest value of N part for the most central bin 
of events and the smallest value of N par t for the most pe- 
ripheral bin of events. On the other hand, the centrality 
selection from the impact parameter, shown by the solid 



curve with no symbols, has the smallest value of (b) for 
the most central bin of events and the largest value of (b) 
for the most peripheral bin of events. When the number 
of charged particles is used to determine the centrality, 
the correlation between the centrality and N part is weaker 
than in the case of using N part . This is especially true 
when only part of the phase space is included, as using 
N c h within — 1 < rj < 1 for the centrality selection (curve 
with open squares) leads to a narrower range of N part 
than using the total number of charged particles for the 
centrality selection (curve with filled squares). The error 
bars shown for two of the curves in Fig. ^| (the curves 
from using N part or impact parameter for the central- 
ity selection) correspond to the root-mean-square widths 
of the impact parameters in each centrality bin, and we 
see that the width in the case of using N part for deter- 
mining the centrality is be quite large (around 1.2 fm) for 
these centrality bins. We also find that, when N c h within 
— 1 < rj < 1 or the total N c h is used for the centrality se- 
lection, these widths are even larger (between 1.3 and 1.8 
fm), indicating that the correlation between N c h and the 
impact parameter is weaker than that between N part and 
the impact parameter. 
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FIG. 14: Charged particle pseudo-rapidity distributions at six 
centrality bins of - 10 - 20 - 40 - 60 - 80 - 100% when the 
centrality is determined from a) the impact parameter or b) 
N c h within — 1 < rj < 1. Dashed curves represent results when 
the centrality is determined from N par t. 

Different centrality selections also lead to differences 
in both the shape and magnitude of the pseudo-rapidity 
distributions at the same centrality bin. Solid curves in 
Fig. 114b and 114b show the pseudo-rapidity distributions 
of charged particles using, respectively, the impact pa- 
rameter and N c h within — 1 < rj < 1 for determining 
the centrality. The distributions corresponding to using 
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Npart for the centrality selection are also shown by the 
dashed curves for comparison. From Fig. 114b . we find 
that dNch/drj at rj = grows slower with centrality when 
the impact parameter instead of N part is used for the cen- 
trality selection. However, it is interesting to see, from 
comparing Fig. 115b with Fig. |SJ that the centrality de- 
pendences of (dN c h/dr}) /Np art for the two cases are quite 
similar. This is due to the compensating effect of slower 
growth of N part with centrality when the impact parame- 
ter is used to select the centrality (see Fig. H^|) . We note 
that in relativistic heavy ion collisions where the mul- 
tiplicity is much higher, these different methods for the 
centrality selection lead to quite similar results, contrary 
to those shown in Figs. 1131 and ^]for d+Au collisions. 
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FIG. 15: Centrality dependence of charged particle multiplic- 
ity with the centrality determined from a) the impact param- 
eter or b) N c h within — 1 < r\ < 1. 



Fig. 115b shows the centrality dependence of charged 
particle multiplicity per participant at different pseudo- 
rapidities when N c h within — 1 < r\ < 1 is used for the 
centrality selection. Comparing with Fig. we find that 
the centrality dependence at r\ = —4 or rj = +4 is similar 
to the case where N par t determines the centrality, but the 
centrality dependence at 77 = —2, or +2 is totally dif- 
ferent. For example, while in Fig. [5] (and in Fig. 115b) the 
curve of (dN c h / drj) /N part at 77 = shows a small decrease 
with centrality, the corresponding curve in Fig. 115b in- 
creases significantly with centrality. We note that these 
differences from Fig. [S] exist even when the total num- 
ber of N c h is used for selecting the centrality (instead 
of N c h within — 1 < rj < 1). Part of this large differ- 
ence in the centrality dependence is due to the stronger 
increase of dN c h/dr](r] = 0) with centrality in the case 
when N c h around mid-pscudo-rapidity is used to deter- 
mine the centrality, as shown in Fig. 114b . Also, since 
the average value of N par t grows slower with centrality 
in this case (see Fig. II 3D . the difference is further en- 
hanced after dividing dNch/dr] by N pa rt- Figs. ITU and 
l*Pol thus show that the criteria of centrality selection can 
introduce appreciable differences in the centrality depen- 
dence of global observables in deuteron-gold collisions, 
and corresponding care must be taken when comparing 
theoretical results with experimental data on centrality 
dependences. 
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